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Abstract 

We study the effects of chiral logs on the heavy — > light pseudoscalar meson 
transition form factors by using standard and quenched chiral perturbation theory 
combined with the static heavy quark limit. The resulting expressions are used to 
indicate the size of uncertainties due to the use of the quenched approximation in 
the current lattice studies. They may also be used to assess the size of systematic 
uncertainties induced by missing chiral log terms in extrapolating toward the physical 
pion mass. We also provide the coefficient multiplying the quenched chiral log, which 
may be useful if the quenched lattice studies are performed with very light mesons. 



PACS:12.39.Fe (Chiral Lagrangians), 12.39. Hg (Heavy quark effective theory), 13.20.-v (Leptonic 
and semileptonic decays of mesons), 11. 15. Ha (Lattice gauge theory). 



1 Introduction 



Over the past decade, a considerable amount of effort has been put into studying the 
nonperturbative QCD dynamics of the heavy — > light decays. The main target was (and still 
is) the extraction of the CKM matrix element \V u b\. The prerequisite for its determination 
from the exclusive B — > -nlv decay mode is a precise knowledge of the relevant form 
factors. Accurate information on the form factors is crucial also when studying the impact 
of physics beyond the standard model (SM) on the exclusive b — > s£ + £~ modes. 

The fact that the kinematically accessible region of momentum transfers is very large 
(e.g., for B — > sxtv it is < q 2 < 26 A GeV 2 ) makes QCD-based calculations of form 
factors ever more difficult. The physical pictures emerging at the two extremities of the q 2 
region are quite different and effective field theory approaches, based on the appropriate 
use of the heavy quark expansion, have been developed to simplify the description of 
these processes. The heavy quark effective theory (HQET), which is applicable for recoil 
momenta close to zero (q 2 — > g^), 1 provides us with valuable scaling laws for the form 
factors [1]. In the region of large recoils (q 2 — > 0), instead, the large energy effective theory 
and its descendants help resolve the heavy quark dependence of the form factors [2,3]. 
Although these conceptual steps forward are highly beneficial for a better understanding 
of the underlying dynamics, a model-independent description (calculation) of the form 
factors in the entire physical region is still missing. 

Among the QCD-based approaches employed to compute the heavy — > light decay form 
factors, the following two stand apart. 

® Light cone QCD sum rules (LCSR). This analytic approach contains the least number 
of assumptions and has the correct heavy quark mass scaling properties. The range 
of applicability is, however, limited to low g 2 's [4]. 

® Lattice QCD. This method allows us to solve the nonperturbative QCD effects nu- 
merically. Because of the current insufficient computing power, the B — > n transition 
is reached either (a) by extrapolating the directly computed form factors from the re- 
gion around the charm to the 6-quark mass [5, 6], or (b) by using a latticized effective 
theory of the heavy quark, such as NRQCD [7] (see also ref. [8]), or (c) by reinter- 
preting lattice QCD in terms of NRQCD when the heavy quark mass becomes larger 
than the lattice UV cutoff [9]. All these strategies share one feature: the accessible 
form factors are restricted to the region of small recoils. 

It is fair to say that the LCSR and lattice QCD are complementary to each other; it is 
important to use them both in order to check their consistency and from their comparison 
perhaps learn more about the underlying nonperturbative QCD dynamics. 

Since the lattice studies are expected to provide us with the most accurate predictions 
about the shapes and absolute values of the form factors, it is important to have good con- 
trol over the various assumptions that are currently used in lattice simulation and the data 
analysis. Two sources of systematic uncertainty have so far been ignored: the quenched 

1 Zero recoil refers to the recoil of the daughter meson in the rest frame of the decaying one. In B — > wiis, 
it means that the pion is soft. 
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approximation and possible deviations from the linear or quadratic chiral extrapolation 
forms. 

All the available lattice results for B — > irlv decay form factors are obtained from 
simulations in the quenched approximation, where the sea quark loop effects in the QCD 
vacuum fluctuations are neglected (tif = 0). To get an idea about the systematic error 
induced by the quenched approximation, one can confront the expressions for the form 
factors derived in the standard and in quenched chiral perturbation theory (ChPT and 
QChPT, respectively). Such expressions, at leading order in the heavy quark expansion 
and next-to-leading order (NLO) in the chiral expansion, are provided in the present paper. 
These expressions are also useful in assessing the systematic uncertainties due to chiral 
extrapolations. Current lattice studies deal with light mesons of masses > 450 MeV. The 
physical pion mass is reached through a linear or quadratic extrapolation in the light 
quark mass. Although it is not clear for which light quark masses one begins to probe 
the subtleties of the chiral expansion, it is beyond reasonable doubt that, very close to 
the chiral limit, chiral log terms of the form m\ log(m^) may modify the result of the 
extrapolation. The coefficients multiplying the chiral logs are predicted by (Q)ChPT and 
will be presented in this paper. Finally, in the case of B — > K decay the standard lattice 
strategy is to consider the kaon as a meson consisting of degenerate quark masses. The 
impact of nondegeneracy in the quenched approximation may also be addressed by using 
the QChPT expressions for the form factors, as we shall see later on. 

It is important to stress a difficulty in getting reliable numerical estimates from this 
approach. As we just mentioned, it is not clear for what value of the light meson mass 
the chiral logs become relevant. That ambiguity is important since one extrapolates from 
the heavier light masses, for which chiral logs have not been observed. Another obstacle is 
the multitude of low energy constants that appear in the Lagrangian and in the transition 
operators in both quenched and unquenched ChPT. The values of some of those constants 
are unknown or simply guessed. Furthermore, as we shall see, the appearance of the chirally 
divergent quenched chiral logs obliges us to compare the full and quenched expressions for 
not-so-light mesons, for which the ChPT is less predictive. For these three reasons, the 
numerical results inferred from this approach should always be taken with a grain of salt. 
In other words, rather than true estimates of the quenching errors, our numerical results 
should be considered as a mere indicator of the size of those errors. 

The rest of the paper is organized as follows. In Sec. 2 we recall the basic elements of 
(Q)ChPT and its combination with the leading order HQET Lagrangian; in addition to 
standard definitions of the form factors, we will introduce some that are more convenient for 
our purposes; in Sec. 3 we present the one-loop (standard and quenched) chiral corrections 
for our form factors; in Sec. 4 we discuss the values of the low energy constants that we 
chose for the numerical analysis, which we present in sees. 5 and 6; we conclude in Sec. 7. 
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2 Setting the scene 



In this section we recall some basic features of QChPT, as it was developed in the pa- 
pers by Sharpe [10] and by Bernard and Golterman [11]. 2 Although QChPT resembles 
standard ChPT in many aspects there are important differences. The main one is the pres- 
ence of the light i]' state, which in quenched QCD (QQCD) does not decouple from the 
octet of pseudoscalar mesons. As a result the "pion" propagator exhibits not only a pole 
structure but also a double-pole one, which is the source of the pathology of the quenched 
approximation. 



2.1 Quenched chiral Lagrangian 

In QQCD, in addition to the quarks q a (a = 1,2,3), one also has the bosonic "ghost" 
quarks q a , of spin | and with identical mass m q ~ a /m qa = 1. Their role is to cancel the 
contributions of the closed quark loops, i.e., they provide quenching. If one assumes that 
the symmetry breaking pattern of QQCD is similar to that of the full QCD, i.e. the 
graded SU(3\3) L <8> SU(3\3) R spontaneously breaks down to SU(3\3) V , the following chiral 
Lagrangian can be written 

flight = ^ str (c^Ec^Et) + ^ str(.ME + ME*) + a d^ d^ - m 2 $ 2 + £ 4 , (1) 
o 2 

where we adopt the convention that / ~ 130 MeV, and the notation 

E = expf 2 4V *=(*f\. (2) 



The following comments are in order. 

- In addition to the standard (qq) Goldstone bosons (n,K,rj) and the vj meson, 3 , 
organized in the 3x3 matrix 



^>°+vV?+vV/ - + K+ \ 



(3) 



the ghost-ghost (qq) bosons (0), as well as the pseudoscalar fermions qq and qq 
(x), also appear in eq. (2). 

The global symmetry SU(3\3) L (g> SU(3\3) R is graded and in eq. (1), instead of the 
familiar trace, one deals with the supertrace str($) = tr(0) — tr(0). 

As already stressed, rj' does not decouple from the light pseudoscalar octet. Its effect 
is included in two terms of the Lagrangian (1), each of them multiplied by a new low 
energy constant, namely, a and m . Note that $ = str($)/v^6 = (rj' — r}')/y/2- 



2 For an elegant alternative way to introduce partially quenched ChPT, see ref. [12]. 

3 Wc will neglect the mixing of rj and ?/ states, as it is irrelevant for the discussion that follows. 
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- The quark-ghost mass matrix is diagonal M = di&g(m u ,m d ,m s ,m u ,m d ,m s ). After 
diagonalizing the mass term in eq. (1), one gets 

ml = 4:fx m q , m\ = 2/i (m q + m s ) , m 2 = (m q + 2m s ) , (4) 

verifying the familiar Gell-Mann-Okubo relation 4m 2 K — m\ — 3m 2 = 0. Notice that 
we neglect the isospin symmetry breaking, i.e., we set m u = m d = m q . 

- In eq. (1) £4 stands for the terms of 0(p 4 ), of which we write only those that are 
relevant to the heavy-to-light form factors, namely, 

£ 4 = str (d^E*) str (ME* + EM*) 

+ L 5 str [d^d^E (ME* + EM*)] +•••}. 

L 4 and L 5 generate the mass correction to the decay constants (and to the wave 
function renormalization constant) [13]. 

It is straightforward to verify that after setting str — > tr, $ — > 0, rj' — > 0, eq. (1) leads to 
the standard (full QCD) chiral Lagrangian. 4 

2.2 Incorporating the heavy quark symmetry 

QChPT has been combined with the leading order HQET in the work by Booth [15] and 
by Sharpe and Zhang [16]. They applied the approach to compute the heavy-light decay 
constants, the B°-B° mixing parameter and the Isgur-Wise function 5 . 

To devise a Lagrangian for the heavy-light mesons, it is necessary to include the heavy 
quark spin symmetry. This is achieved by combining the pseudoscalar (P a ) and vector 
(P* a ) heavy-light mesons in one field: 

H a (v) = ±±£ [P; » 7 , - P»7 5 ] , (6) 

where (l+^)/2 projects out the particle component of the heavy quark only. The conjugate 
field is defined as H a (v) = 7 iJ](t>)7 , whereas the covariant derivative and the axial field 
have the following forms: 

DlH h = d»H a - H b Yl = d»H a - H^d.t + , 

A? =^K t ^-^ t U- (7) 

In the above equation, a and b run over the light quark flavors and £ = \/E. With these 
ingredients in hand, we now write the quenched chiral Lagrangian for the heavy-light 
mesons as 

- str a Tr[H a ivD ba H b ] + g str a Tr[H a H b ^A^ 5 ] 

4 For a review of the standard ChPT see one of the references listed in [14]. 
5 For a recent result on the Isgur-Wise function in partially QChPT, see ref. [17]. 
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+$f str a Ti[H a H aWj ] str(A^) + £ 3 , 



(8) 



where g (g') is the coupling of the heavy meson doublet to the Goldstone boson (to rj' or 
fj'). A term with g' is thus an artifact of the quenched theory. The higher order terms in 
the expansion in v-p ~ 0(p), 0{p 2 ), and in m q ~ 0(p 2 ), denoted as £ 3 in eq. (8), have the 
following form [15] 

C 3 = 2Ai str a Tr[H a H b ] (M+) ba + h str a Tr[H a ivD bc H b ] (M + ) ca 

+k 2 str a Tr[H a ivD ba H b ] str c (M + ) cc + . . . , (9) 

with M+ = (l/2)(gM& + £M£). We displayed only the terms that contribute to the 
heavy-to-light form factors. In the above equations, "Tr" stands for the trace over Dirac 
indices, whereas "str" is the supertrace over the light flavor indices. 

As in the previous subsection, one can easily verify that after replacing str — > tr, $ — > 0, 
rj' — > 0, and g' — > 0, one recovers the standard chiral Lagrangian for heavy-light mesons 
(for recent reviews see ref. [18], and for the original papers see ref. [19]). 



2.3 Form factors 

A frequently encountered decomposition of the matrix elements relevant to the leptonic, 
the semileptonic and the penguin-induced hadronic matrix elements is 

{0\qi^l5b\B(p B )} = if B p Bf , , 

(P(p)\qi,b\B(p B )) = [{pB+p), - q ™*~ m2p > ) F + (q 2 ) + ^^ W) , 

(P(p)\q^uq u b\B(p B )) = i ( q \p B+ p)^-( m l- m l) q \ , (10) 

V / m B + m P 

where q = d or s, \P{p)) is the light pseudoscalar meson state (pion or kaon), and q u = 

(pb-pT- 

We will be working in the static limit, m B — > oo. The eigenstates of the QCD and 
HQET Lagrangians are related as 

lim p — \B(pb))qcd = |5(v)) H qet • (11) 

m s -+oo ^Jm B 

In the static limit it is more convenient to use definitions in which the form factors are 
independent of the heavy meson mass, namely, 

(0 1 97^,7561; I -B(u)) hqet = ifvp , 

(P(p)\q~l»b v \B(v)) mET = [pp- (v-p)v^]f p (v-p) + v^f v (v-p) , (12) 

where the field b v does not depend on the heavy quark mass [20]. The form factors f PjV are 
functions of the variable 

m% + m 2 P -q 2 

V ' P = 2^ ' (13) 

6 



which in the heavy meson rest frame is the energy of the light meson Ep. The relation 
between the quantities defined in eqs. (10) and (12) is obtained by matching QCD to HQET 
at the scale // ~ m b [21]: 



= C ll (m b ),Jrn B ~[f p (v-p) + 0(l/m B )] , 
(m B + E P )F + (q") - (m B - E P ) ml ~f p [F + (g 2 ) - F (g 2 )] |^ 

= C J0 (m b )y/W^f v (E P ) + 0(l/m B ) , 

2m B 2 

-F T (q ) = C lll0 {m b )f p {v-p),Jrn B ~ + 0(l/m B ) ■ 



(14) 



In the following, we set the matching constants C 7i to their tree level value (C 7 = 1). By 
neglecting the subleading terms in the heavy quark expansion, we have 6 

F °(<1%^ =-±=Mv.p), F + (q% =F T (q\ =^f P (vp), (15) 

which exhibit the usual heavy mass scaling laws for the semileptonic form factors [1]. The 
equality F T (q 2 ) = F + (q 2 ) arises from 

(P(p)|^°(l + 75)6«|5(U)> H QET = - 75)6«|5(U)> H QET , (16) 

which is easily obtained by using the equation of motion 7 &„ — b v . In the following our 
main concern will be the evaluation of the long-distance (chiral) corrections to the form 
factors f v ,p(v-p). 

2.4 Heavy- light current 

Having in mind eq. (16), we only need to consider the (V — A) Dirac structure. In the 
static heavy quark limit and at next-to-leading order in the chiral expansion, the bosonized 
heavy — > light current reads [15] 

r = q a r(i - 7 5 )q - y - ^mli 1 + vm%) 

+ Tr[ 7 ^(l - TBWttiM+U (17) 

+ y x 2 Tr[ 7 ^l - l5 )H b ]Zl str c (M + ) cc . 

When bosonizing the current J M , one could envisage an arbitrary function Vl(^o) m the 
first term on the right of eq. (17), expanded in a Taylor series. Only the terms linear in 

Corrections at 0{\/mb) have been discussed in [22-24]. 
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Figure 1: The point (1) and the pole (2) tree level Feynman diagrams contributing to heavy — ► light 
transition form factors. The box denotes the weak current insertion. 



$o are relevant to our purpose and Vj,(0) can be normalized to 1 [16]. The appearance of 
V[(0) is yet another artifact of the quenched theory. The phase of the heavy meson can be 
chosen in such a way that V[(0) is completely imaginary, whereas the constants a, xi, and 
x 2 are real. Notice that, at leading order in the chiral expansion, the constant a is simply 
the heavy-light meson decay constant in the static limit (mg — > oo), i.e., / in eq. (12). 

3 Chiral Corrections to f p (v-p) and f v (v-p) 

In this section we explain the main steps undertaken to compute the chiral logarithmic 
corrections for /, f p (v-p), and f v (v-p). Our results for the decay constant / agree with those 
of refs. [15, 16]. They will be used in Sec. 5.3 to construct the ratios that are independent 
of counterterms. In ref. [25] ChPT was applied to compute the heavy-to-light form factors. 
We repeat that calculation and extend it to the quenched case. 

3.1 fp rec (v-p) and fj iec {v-p) 

We start the discussion by writing the tree level expressions 

f P (vp) = J v . p 9 +A * ' f*( v 'P) = J ' t ' = a ' ( 18 ) 

The point and the pole diagrams, depicted in fig. 1, describe /„ and f p , respectively. The 
index V in A* = m^* - wb, labels the light quark flavor. When necessary, we will use 
Pij to denote the light pseudoscalar mesons with valence quark content q^. 

Although the heavy quark spin symmetry suggests A* — > 0, we will keep this term 
finite because it provides the pole at m 2 B * to the form factor f p , or equivalently to the form 
factors F +jT . 7 The form factor /„ (or F ), on the other hand, does not depend on (v-p) at 
the tree level. 

3.2 f° ne - lo °v(vp) and /° nc - lo °P(^j9) 

Neglecting the crosses, all the graphs shown in fig. 2 describe the one-loop chiral contribu- 
tions to the form factors f v , P (v-p) in the full ChPT. In the quenched theory, in contrast, 

7 The pole dominance is easily seen if one rewrites the denominator of f p (v-p) as v-p + A* = 
(ms»/2) (l — q 2 /m%,) 7 where the corrections (mj. — ms)/ms and m 2 P /m 2 B are neglected. 
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Figure 2: The one- loop contributions to the B — > 7r transition. Double (single) lines denote the heavy 
(light) meson, while the weak current insertion is depicted by the empty box. Crosses stand for the m 
hairpin vertex (1). Each graph represents two Feynman diagrams in quenched ChPT: (a) the diagram 
without the hairpin vertex (cross) and (b) the same diagram with a cross. In full ChPT only the diagrams 
without the cross are present. The amplitudes corresponding to the diagrams are listed in appendix B. 



graphs both with and without crosses appear. The cross denotes the so-called hairpin 
vertex, i.e., the dipole term in the "pion" propagator: 

1 ( 1 M _ a °p 2 - < ( 1 i \ m 

p 2 -m 2 P \0 -I J (p2 _ \ 1 1 J ■ 

In the computation of the loop integrals, the naive dimensional regularisation has been 
used, with the renormalization prescription MS + 1, where A = 2/e — 7 + ln(47r) + 1 is 
subtracted [13]. We neglect the isospin symmetry breaking {m u = m d = m q ), as well as 
the mass differences between B, B s , B*, B* meson states, whenever they appear in the 
loop. The resulting expressions can be written as 

Quenched ChPT 



f^ p Hvp) = i 



a 



9 



f v-p + A* 



1 + bfp^ Pl3 + ^-rn-j + kwh - AL^-irm + mj ) 



a 



1 + Stf*^ + hf + ^ )mj - 4L 5 ^( mi + rrij) 



P 



fj = a 1 + Sfj + [ y + X! ) m 3 



(20) 
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Full ChPT 



a g 
f v-p + A* 



1 + 5f; V 3 + — m j + K i m i - 4L 5 -^(mi + rrij) 



k 2 _ 

2 



- ( — - 8L 4 — + k 2 ) (m u + m d + m s ) 



a 



1 + Stf^ + (| + xx) m 3 - 4L 5 ^K + m 3 ) 



. k 2 T 4/i ' 



1 + + ( y + x i J m i + ( y + x 2 ) (mu + m d + m s ) 



(21) 



We separated the chiral loop corrections (Sf) from those involving the counter-terms. The 
loop corrections to the form factors are written as 



(22) 



where the sum runs over all diagrams shown in fig. 2, and SZb^p^) encodes the chiral 
loop contributions arising from the wave function renormalization of the heavy-light and 
light-light meson respectively. The explicit expressions are lengthy and are collected in 
appendixes B.l (for the quenched) and B.2 (for the full theory). As can be seen from 
eqs. (20) and (21), no dependence on (v-p) arises from the counterterms. The modification 
of the tree-level (v-p) dependence of the form factors f v>p is entirely due to the chiral loop 
corrections. It is worth noticing that in the quenched approximation the form factor fy^ n 
is completely independent of (v-p). This is so because the contribution from diagram 4 in 
fig. 2 vanishes in the isospin limit m u = ma- 

An important feature emerging from this calculation is that the quenched form factors 
f VtP and quenched decay constant / suffer from the common quenched pathology, that is, 
from the presence of the quenched chiral logs oc ml \og(m 2 P / fj 2 ). Such terms are divergent 
in the limit mp — > 0, suggesting that in the quenched approximation the chiral limit for 
any of F?~^% or is not defined. This feature is also present in the light meson sector, 
e.g., for the chiral condensate, the light meson decay constant consisting of nondegenerate 
quark flavors, etc. [26,27]. 



4 Choice of Parameters 

For numerical analysis of the expressions given in eqs. (20) and (21) we need to make a 
specific choice of quite a number of low energy constants. To do so we will mainly rely on 
the existing lattice data. In table 1 we collect the parameters whose (range of) values we 
were able to fix. In the following we briefly discuss each of those values. 
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Coupling Full theory Quenched theory 



a (GeV 3 / 2 ) 


0.56 ±0.04 


0.48 ±0.03 


9 


0.50 ±0.10 


0.56 ±0.12 


f (GeV) 


0.130 


0.124 ±0.004 


Ho {(jev ) 


1 i/i xn on 
1.14 It U.ZU 


1. lo It U.U4 


L A (xl0~ 3 ) 


-(0.5 ±0.5) 


(0.0 ±0.5) 


L 5 (xlO- 3 ) 


0.8 ±0.5 


0.8 ±0.2 


m (GeV) 




0.64 ±0.06 


9' (GeV) 




-0.6 to 0.6 



Table 1: Low energy constants whose determination is discussed in the text. 

o oc. This constant is equal to the heavy-light decay constant in the static limit of QCD. 
Its value can be obtained from fo\^riD and/or /s^/ms, which are then extrapolated 
(in inverse heavy-light meson mass) to the infinite mass limit as 

f H V^ = a(l + A/m H + B/m 2 H + ...) . (23) 

From lattice QCD and QCD sum rules, it is known that the slope A is large and 
negative, whereas the value of the curvature B is small. Therefore, to a good accuracy, 
one can set B = and neglect higher terms in l/m H . In this way we obtain the 
following: 

— From extensive quenched lattice simulations by the MILC collaboration, one 
can deduce a qucnch = 0.45(5) GeV 3//2 [28]. In the same paper, they also present 
the results of their unquenched simulations, from which we extract o/ u11 = 
0.53(7) GeV 3 / 2 . These numbers agree quite well with the values obtained 
by CP-PACS Collaboration, namely, a ^ cnch = 0.50(4) GeV 3/2 , and a fu11 = 
0.57(6) GeV 3/2 [29]. Notice that both references use the same treatment of 
the heavy quark on the lattice (the so-called Fermilab formalism). 

- The UKQCD [30] and APE [31] Collaborations compute heavy-light meson de- 
cay constants using the fully relativistic lattice QCD, but only for the mesons 
of masses run £ (1.8,2.7) GeV. From a linear fit of the form (23), from their 
quenched data UKQCD obtain 8 a qucnch = 0.49(4) GeV 3/2 , in agreement with 
the previous result by APE, a <i ucnch = 0.48(5) GeV 3/2 . 

8 We thank G. Lacagnina for communicating this result to us. 
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— Recent results obtained by using the QCD sum rules agree quite well with the 
above unquenched values. From the compilation of the QCD sum rule estimates 
in ref. [32], one has a fu11 = 0.58(9) GeV 3/2 (for the most recent result see [33]). 

o g: This constant is related to the phenomenological coupling gD*Dn as 

2y/m D m D * 

9d*Dit = 7 9 , ( 24 ) 

where one can also set mo* = tub because the above relation is defined only in the 
static limit, in which the heavy quark spin symmetry is exact. From the experi- 
mentally measured total width of the D* + meson, one gets ^ u11 = 0.59(7) [34]. The 
subscript "c" warns us that the value is obtained in the charm quark mass sector. On 
the lattice, that value has been computed recently in the quenched approximation, 
leading to g^ uench = 0.66(9) [35]. To get the value of g one needs to extrapolate to 
the infinite heavy quark mass limit, i.e., Qq = g + j/itih + ■ ■ ■ The lattice data of 
ref. [35] suggest that the slope 7 is negligible, and thus g/g c ~ 1. On the other hand, 
the LCSR calculation suggests that the same slope is significant and positive [36]. 
By neglecting the terms of 0{l/m 2 H ) and higher, ref. [36] leads to g/g c ~ 0.7. We 
will take the average of the two (lattice and LCSR), that is, g/g c « 0.85, and add 
the difference in the error bars of both quenched and unquenched values. 9 

o /: To get the chiral coupling constant we will use f n = 132 MeV and fx = 160 MeV, 
and the fact that m s /m q = 24.4 ± 1.5 [38]. After linearly extrapolating to the chiral 
limit (m q — > 0), we get / = 130 MeV. The recent extensive quenched lattice study 
with Wilson fermions gives / = 119 ± 7 MeV [39], whereas the one with staggered 
fermions gave / = 125 ± 9 MeV [40]. The latter result was also obtained on larger 
lattices with domain wall fermions, namely, / = 125 ± 7 MeV [41]. As a weighted 
average, we will take / = 124 ± 4 MeV. 

o fi : From the Gell-Mann-Oakes-Renner formula [42], it is easy to identify fi^ = 

— (qq)/f%- Its value (in the full theory) can be fixed by using (qq) MS (2 GeV) = 

— [267(16) MeV] 3 [43]. The quenched value is extracted from the lattice data. Recent 
results in the MS scheme (at the renormalization scale 2 GeV) are // = 1.13(4) GeV [39 
and fi = 1.10(11) GeV [44]. With these numbers and by using / = 124±4 MeV, the 
corresponding numerical values for the chiral condensate in the quenched approxima- 
tion are (qq) m (2 GeV) = -[259(6) MeV] 3 and (gg) m (2 GeV) = -[257 (10) MeV] 3 . 

° £4,5: These two couplings have been extensively discussed in the literature [14]. Their 
values in the full theory, at the scale /i ~ 1 GeV, are given in table 1. In quenched 
QCD, the coupling L 5 was recently computed in ref. [45], with the result a$ = 0.99(26) 
and L 5 = a 5 /(1287r 2 ) = 0.8(2) x 10 -3 . On the other hand, the quenched estimate of 



9 We are aware of the result of ref. [37], g — 0.42(9), where this coupling was computed for the first 
time on the lattice, in static HQET. However, in view of the insignificant statistics, very coarse lattice, 
and only two light quark masses, the value g = 0.42(9) should be considered as exploratory. An improved 
calculation of g, along the lines described in ref. [37] would be most welcome, though. 
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L 4 is not available. On the basis of the large N c expansion we only know that L 4 is 
smaller than L 5 (see ref. [13]). We will take it to be zero and vary by ±0.5 x 1CT 3 . 

o m : This mass characterizes the magnitude of the hairpin insertion (crosses in fig. 2). 
It enters in the coefficient of the quenched chiral log, which, in the literature, is often 
referred to as 5 = m\j \24ti 2 f 2 ). The precise value of 5 is unknown at present, mainly 
because in realistic lattice studies it is very difficult (if possible at all) to resolve 
the finite lattice volume effects from the quenched chiral logs. This is why it is still 
not completely clear whether the observed deviations from a linear dependence of 
the pseudoscalar meson squared {nip) on the light quark masses (m q ) is due to the 
presence of the chiral logs, or if it is an artifact of the lattice geometry (see, e.g., 
ref. [46]). With this remark in mind, we now quote recent values for the parameter 
5, as obtained from lattice studies of the dependence of mp on m q : 

- with Wilson fermions, 5 = 0.10(2) [39]; 

- with modified Wilson fermions, 5 = 0.073(20) [47]; 10 

- with domain wall fermions, 5 = 0.029(7) [48]; 

- with overlap fermions, 5 = (0.0-0.2) [49] and 5 = (0.2-0.3) [50]. 

For the pre-1996 results see the review in [27]. The value of 5 obtained by the CP- 
PACS Collaboration stands out because they made an extensive high statistics study 
on large lattice volumes and extrapolated to the continuum [39]. A worrisome aspect 
of this value, however, is that the light pseudoscalar mesons that they access directly 
lie in the range m P G (300,750) MeV, for which the significance of the chiral logs 
may be questionable. Keeping this comment in mind and using their 5 = 0.10(2), we 
get m = 0.64(6) GeV. 

o g'\ The value of this constant can be obtained from Qbb*t\' and <7ddv the couplings 
of the lowest lying heavy-light meson doublet with the light quenched rj' state. Such 
a lattice study has not been made so far. To get a rough estimate we may rely on 
the large N c limit from which one expects \g'\ < g. 

As for the other parameters appearing in the quenched expressions, in the numerical anal- 
ysis, we will first set them to zero and then vary their values in the ranges suggested by 
the large N c expansion. More specifically, we take \V[\ < 0.5, |xi|, |fci| < 32fi L 5 / f 2 , and 
|^2| = I&2I = 0. The parameter a will be varied in the range | ct 1 < 0.1. It is important 
to stress that the effect of the variation of these latter parameters is completely negligible 
for our numerical estimates. 

Note that the low energy constants in table 1 are extracted from lattice data at the 
tree level of ChPT (except for the counterterms) . 11 

10 Thc final result of ref. [47], S = 0.065(13), is obtained by combining various ways of extracting this 
quantity from the lattice data. For an easier comparison with other groups, we quote the result given in 
eq. (29) of ref. [47], which is obviously fully consistent with their final value. 

11 The standard procedure would involve "undressing" the chiral loops from the measured low energy 
constants. This procedure is not applicable when using the lattice data because direct lattice computa- 
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Figure 3: The ratios Q P , v (v-p) are defined in eq. (25). The plots refer to the mesons consisting of two 
degenerate (valence) quarks with mass m q — r m s and r = 0.25. The counter-term coefficients fci.2, >fi,2, 
as well as |V£(0)| are neglected, while \i = f GeV ~ A x . Values of other parameters are given in table 1. 

5 Quenching Errors 

In this section we will use the expressions for the form factors (20) and (21), the parameters 
listed in table 1, to get an estimate of the quenching errors. We reiterate that the results of 
this section refer to the zeroth order in the l/rnn expansion and to the leading logarithmic 
chiral corrections, with the specific choice of parameters discussed in the previous section. 
We also stress that in all the following discussion we will keep the strange quark mass 
fixed to its quenched value of mf s (2 GeV) = 105 MeV, and m{ ul1 j m f ench ~ 0.85 [51]. To 
examine the quenching errors we will study the following ratios 



The parameter that makes the strongest impact on the results for Q PtV is g' . As stated 
above, its value is expected to lie in the range —g<g'<g. Since even its sign is not 
known we will distinguish between the two "extreme" cases g' = —g and g' = +g. 

5.1 B — > it transition 

To examine the quenching errors for B — > 7r decay we cannot set the pion mass to its 
physical value, because the spurious quenched chiral logs would become dominant. How- 

tions are made at m n > 400 MeV, and the physical results are obtained through a linear or quadratic 
extrapolation in — > (rn^) phvs . These extrapolations do not include the chiral logs, so "undressing" the 
chiral loops would lead to unrealistic estimates of the low energy constants. For this reason, we adopted 
the tree-level approximation in extracting the low energy constants. 




flf (v-p) - f*? nch (vp) 



(25) 
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ever, we cannot go far away from the chiral limit either, because the (Q)ChPT approach 
then becomes inadequate. To those two competing requirements, we should add a third: a 
desire to be sufficiently close to the region of quark masses (i.e., pseudoscalar meson masses 
squared) probed by the current quenched lattice studies [5-9] . It is not clear whether or not 
a mass interval for which all of the above requirements are satisfied exists. We will assume 
that they are satisfied for trip ps 330 MeV. From eq. (4), a "pion" of mass mp ps 330 MeV 
is composed of two quarks of mass corresponding to r = m q /m s ~ 0.25, with respect to 
the strange quark mass m s . In fig. 3 we plot the ratios Q P)V (v-p,r = 0.25). 

We notice that, regardless of the value chosen for g', the quenching errors in the form 
factor f p (v-p) are not excessively large for most of the v-p where ChPT is applicable. This 
is important since this form factor is the only one entering the B — > -niv decay rate, from 
which we hope to be able to extract the value for 

By specifying g' = +g, we observe that Q p (v-p) > 0, for any (v-p)/ p < 0.7 and r > 0.2. 
In other words, quenched values for f p are smaller than unquenched ones. If we take 
g' = —g instead, the ratio Q p (v-p) has a zero. The point of zero quenching errors in 
Fig. 3 is found at (v-p/p) = 0.47(3) for r = 0.25 and the values of parameters given in 
table 1. Figure 4 shows the curve of vanishing quenching errors in the v-p — r plane, i.e., 
Q p (v-p, r) — ± 10%. This curve really exists for g' = —g, while for g' = +g only the part 
corresponding to Q p (v-p, r) = +10% can be reached for r < 0.4. In summary, it is possible 
to find combinations of the pion mass and the pion recoil energy such that the quenching 
errors in the dominant form factor f p (v-p) are kept under the 10% level. 

^From fig. 3 we also see that there exists a point (v-p/ p) ~ 0.3 such that the ratio Q p 
is independent of the value of g'. At that point, we get 

Q p (0.3p, r ~ 0.25) ~ 20% . (26) 

This result is a useful estimate of the quenching error for realistic values of (v-p) and r, 
which are currently probed on the lattice. For the reader's convenience, we have fitted the 
points corresponding to (df p /dg') = to a polynomial in r e [0.1, 1] and obtained 

— = 0.17(1) + 0.8(l)r - 1.3(2)r 2 + 0.58(2)r 3 . (27) 
p 

We also checked that for 0.1 < r < 0.35, with r and (v-p) satisfying eq. (27), the quenching 
errors are kept under the 25% level. This concludes our discussion of the form factor f p . 

The situation with the form factor f v (v-p) is much worse. From fig. 3, we see that 
the quenching errors are in the range 30 — 60%, and drop below that level only for larger 
recoils for which the present approach is not appropriate. The quenching errors remain 
large when varying r in the range r G [0.1, 1]. A somewhat less pessimistic situation is 
present at zero recoil (v-p = 0) where, contrary to the case of f p , the form factor f v can be 
extracted from the lattice data. For m P = 330 MeV, at zero recoil we get 

f 8% (g' = +g), 
Q v (0,r = 0.25) ~ { 

{ 23% (g> = -g) . 

Once again, we warn the reader that the numerical estimates made in this section are 
obtained for the set of low energy constants specified in table 1. 
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Figure 4: The thick dashed line corresponds to the zero quenching errors [Q p (v-p,r) = 0] in the case 
when <?' = —g (such a line is not accessible for g' = +g). The shaded region corresponds to the variation 
of Q p (v-p,r) by 10%. Notice that in the case g' = —g, the upper (lower) curve corresponds to Q p (v-p,r) 
equal to —10% (+10%). For the case g' = +g, only the region Q p (v-p,r) < +10% is shown. 
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Figure 5: The ratios Q P , v (v-p) from eq. (28). The plots refer to mesons whose one valence quark is fixed 
to the s-quark mass, while the other has the value m q = rm s . The thick gray vertical line marks r ~ 1/25 
for which the physical kaon mass is reached. As before, the counterterm coefficients fci^, xi,2, as well as 
V[(0) are neglected, while fx = 1 GeV. 



5.2 B — > K transition 

We now turn to the case of a kaon in the final state. In the discussion we shall fix one of 
the valence quarks to the strange quark mass, and vary the other one (jn q = r m s ) in the 
range 1/25 <r< 0.5. In this way the "kaon" mass is varied in the range m p ^ ys < mp < 
600 MeV. For simplicity we consider {v-p)/fi ~ 0.3 (reasonably small recoil), and examine 
the quenching ratios 

g , (0 . 3ft ^^«^. (28) 

U,v (V-p) 

In the quenched expressions for f P)V , given in appendix B, we set for the active quark 
Mf = 4/x m s , while for the spectator one we take Mj = 4/i m s r [see eqs. (46), (47)]. As 
in the previous subsection, in fig. 5 we plot Q^ v for the two "extreme" scenarios, namely, 
g' = +g and g' = —g. We again observe that, regardless of the value for g' , the quenching 
errors on the form factor f v (v-p), > 0.5. The quenching errors Q p are only moderately 
large if g' = —g, and are large for g' = +g. 

Therefore, from this and the preceding subsection, we see that in both channels B — » 7r 
and B — > K our approach suggests that the quenching errors in the form factor f v (v-p) 
[i.e., F (q 2 )] are large > 30%. The quenching errors in f p {v-p) [i.e., F +t T(q 2 )}, on the other 
hand, depend crucially on the value of the low energy constant g'\ they are large for g' > 0, 
and "not so large" for g' < 0. 
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Figure 6: The ratios R PtV defined in eq. (29), measuring the errors induced by the degeneracy in the 
"kaon" quark masses in the quenched calculations. 



Before closing this subsection let us mention that in the quenched lattice studies the 
kaon is usually considered as a composite state of two degenerate quarks. Using mp = 
4fiQm s r, one varies r in order to reach the physical kaon mass rap = m P j^ ys , which occurs for 
r ~ 0.5. One may wonder if the form factors with such a kaon differ from the ones in which 
the quarks in the kaon are nondegenerate, with one of the quarks fixed to the strange mass 
and the other varying toward r ~ 1/25 (i.e., the physical kaon mass). To keep the mass 
of the pseudoscalar meson the same in both situations, we will vary r de g. G [0.6, 0.8] and 
r nd g . = 2 r deg. - 1 £ [0.2,0.6]. As before, we take (v-p/fi) = 0.3 and examine the following 
ratio 



Hv,p(j'deg. 



fv d p 9 \v-p,r ndg ) - f^ p 9 iv-p,r deg . 



rndg. 
Jv,p 



v-p, r ndg ) 



(29) 



ndg. 



=2r, 



deg. 



in the quenched theory. From fig. 6, we again observe the same dichotomy: for g' = —g 
the situation is quite favorable, i.e., the errors due to degeneracy in the quark masses are 
very small, whereas for g' = +g the form factor fp e9 ' is highly overestimated with respect 

to/;*-. 



5.3 Useful ratios 

The double ratio f^"^ K Jb/ fv^ K fs s is ver y gratifying from the ChPT point of view. In 
a double ratio of the form 



J y(Q) 

$i > Pj i 

(Q) 



f 



(Q) 



1 + S fv(Q) Pv ~ S fv(Q) J ' + fifBi (Q) - <*/Bj(Q) , 



(30) 
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where Bj ~ bq~j and P^ ~ q^qj, the dependence on a cancels completely and so do the 
counterterms. Quenching errors in this quantity are thus far more reliably predictable in 
the framework of ChPT than for the separate form factors. This quantity might then prove 
useful in future lattice simulations with both B s — > P and B — > P transitions. 

A similar (partial) cancellation of dependence on counterterms occurs also for the 
quenching error in the ratio f p B ^* K / fp T* 



j.Bj—>Pij D. 

Jp,v(Q) Jp,v 

Jp,v(Q) Jp ' v (31) 

8 

P 



m\ , 



which does not depend on counterterms if L^ ul1 = jj^ ench . In eq. (31) we have neglected 
counterterms suppressed by m\jm 2 K and assumed that 5f is small (5f is independent of 
counterterms by construction). 



6 Chiral extrapolation in B —> n decay 
6.1 Quenched case 

So far the quenched lattice studies of the B — > 7r transition matrix element were confined 
to the region of not very light pseudoscalar meson masses which (from the point of view 
of QChPT) is almost fortunate since one stays away from the region in which the spurious 
quenched chiral logs dominate over the other (physical) contributions. Supposing that one 
manages to push the lattice studies toward ever lighter "pions", the quenched chiral log 
will become a more important effect. From the expressions presented in Sec. 3 and in 
appendix B, for the very light meson mp < (v-p), we obtain 



f P (v-p) 



ag 



f(v-p + A*) 



1 + (cf + <%m 2 p) In (^f ) + eg + c\m P + c p 2 m 2 P + ... 



a 



fv(v-p) = J 



1 + (cf + 4ml) In (^f ) +c v + c v 2 m 2 P + ... 



/i 



(32) 



The quenched chiral log term is proportional to 

m 2 (l + 3g 2 ) 



Si 



6(4tt/)< 



(33) 



and it is a quenched artifact that diverges in the chiral limit. The presence of this term has 
an important consequence, that by approaching ever lighter mp, the form factors fQ uench 
and fv uench increase. This is in contrast to the full (dynamical) theory in which the effect 
is the opposite, i.e., the chiral logs lower the form factors for small masses, as we shall see 
in the next subsection (see fig. 7). The form factor f p (v-p) also picks a term linear in m P in 
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this limit, which is yet another artifact of the quenched approximation. The accompanying 
coefficient reads 



g 2 m\ Air 



(4tt f) 2 3vp ' 

6.2 Unquenched case: An illustration 



(34) 



The unquenched equivalents to the expressions (32) are not straightforward to derive. The 
main obstacle comes from the fact that, instead of one mass in the integrals J 2 (M, v-p) 
and Ji(M,v-p), one now has three masses: M G {m^, m K , m v }. The behavior of those 
integrals depends on the sign of 1 — (v-p)/M (see appendix A). The variation of the pion 
mass entails a change of rrix and m v , which straddle (v-p). This then changes the behavior 
of the integrals J 2 and Ji. For the B — > n transition, we have 



f P (v-p) 



ag 



f (v-p + A* 

l + 3g 2 



1 + 



1 



g 2 ^4:J 1 (m n ,v-p) + 2Ji(m K ,v-p) + ^Ji(m^,f-p)^ 



(4tt/) 2 

, , y 9ml log(m 2 ) + Qm 2 K \og(m 2 K ) + m 2 log(m 2 ; 



dl 



(35) 



fv(v-p) = 



- i 1 

71 + WfT 2 



15-27# 2 2l . 2 . l-3g 2 2 , . 2 . 

■m ff log(mJ + — - — m K \og(m K ) 



12 



1 + 3^ 
12 



m 2 log(m 2 ) + 2I 2 {m. K ,v-p) + I 2 (m K ,v-p) 



+ dg + ... , 



where d is a constant and the ellipses stand for the higher order terms in the rri^ Kr) 
expansion. 



(v-J>) 


^(GeV 1 / 2 ) 


a«(GeV- 3 / 2 ) 


^(GeV" 1 / 2 ) 


4 1} (GeV- 5 / 2 ) 


Reference 


0.55 GeV 


2.5(2)1°:° 


i-i(2)18:§ 


o.9(i)i83 


0.7(1)18;? 


[6] 


0.19 GeV 


0.8(3) 


0.3(3) 


4.8(4) 


3.6(4.5) 


[7] 



Table 2: The parameters describing the linear chiral extrapolation for the B — > 7r form factors 
fl a p(v-p) at fixed (u-p), as indicated in eq. (36). The values are those obtained in the quenched 
lattice studies in refs. [6, 7]. 

To exemplify the impact of the chiral logs, we will now use the existing (quenched) 
lattice results for the B — > n form factors presented in refs. [6] and [7] , in which the chiral 
extrapolation has been made linearly, i.e., 



/£?(«*) = O-P) + a^(vp)m 2 P , 



,(o) 



(36) 
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The parameters civ,p of ref. [7] are obtained by fitting in the region of light pseudoscalar 
mesons that corresponds to mpG (450, 800) MeV. Those of ref. [6] are obtained from a fit 
in mp E (540, 840) MeV. The numerical values are given in table 2. 12 

We will now assume that (a) these results are the same in the full (unquenched) the- 
ory; 13 (b) the form (36) holds true down to a point itlm ~ 250 MeV or Tom k 330 MeV, 
where we smoothly match the lattice and ChPT results. The full ChPT form factor, given 
in appendix B, is used from the matching point itlm down to the physical pion mass. In 
other words, for a fixed value of (v-p), we take 

U,Ml) = 0(m 2 p-m 2 M)f p y(m 2 p) 



+0(m M - m 1 



(dfSf T /dm P \ m%=mli - «S) K - ™m)] • (37) 



In fig. 7 we show the effect for the form factor f p (v-p). We observe that the form factor 
obtained by including the chiral logs in the extrapolation to mp = m\ is smaller than 
the one obtained by extrapolating linearly (dotted lines in the plot). The amount of that 
suppression obviously depends on the value of m 2 M = 4/j,om s rM' the effect of the chiral 
log is smaller for smaller tm- In our example we took tm = 0.15 (rriM — 250 MeV) and 
r M = 0.25 (m M ^ 330 MeV). Using 

f eq.{37) _ f eq.(36) 

«</*.) = - f .,<ir ■ < 38 > 

JP,v 

to measure the error due to chiral logs that were not included in the extrapolation to the 
physical pion mass, we obtain the following results. 

r M = 0.15 

err(f p (vp = 0.19 GeV)) ~ -2% , err(f p (vp = 0.54 GeV)) ~ -5% , 



r M = 0.25 

eri(f p (vp = 0.19 GeV)) ~ -7% , en(f p (vp = 0.54 GeV)) ~ -15% (39) 

The above analysis applied to the form factor f v leads to even smaller uncertainties: 
err(/„) r . M=0 .i5 < 3% and err(/„) rM= o.25 < 6%, for both values of (v-p). 

This exercise is made just to illustrate how one can proceed in order to get an estimate 
of the systematic uncertainties due to the chiral extrapolation. As we saw, the amount of 
estimated uncertainty is highly sensitive to the choice of the point r M , which is why the 
outcome of this exercise should be considered only as a rough estimate. 

It is important to stress again that had we used the quenched expressions (32) to guide 
the chiral extrapolation, the result would stay above the result of the linear extrapolation, 
precisely the opposite to what happens in the full theory (which we show in fig. 7). 

12 We are particularly indebted to Tetsuya Onogi for communicating these results to us. 
13 Since the purpose of the discussion in this section is to illustrate the impact of the chiral logs on the 
result of an extrapolation to the physical pion mass, this assumption should not worry the reader. 
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Figure 7: Illustration of the effect of inclusion of the chiral log in extrapolating to the physical pion from 
the masses accessed from the lattice simulations (marked by the shaded boxes). The full (dashed) curves 
are obtained by including the chiral logs in extrapolation starting from Tm = 0.15 and from 0.25 (gray 
vertical lines) using full (unquenched) ChPT. For comparison we also show the (dotted) line corresponding 
to the linear extrapolation. 



Finally we would like to make a comment on the ratio f v /{f/f) or equivalently F (v- 
p)/(/b//it), which, according to the soft pion theorem, should be equal to 1 at the zero 
recoil point v-p — * 0. In the quenched theory the form factor Fq is independent of v-p and 
its chiral corrections are exactly the same as in fs/fn- Therefore the ratio is indeed 1 for all 
combinations of (r,v-p). In the full (unquenched) theory, in contrast, the chiral corrections 
cancel only at v-p — > and the soft pion theorem is satisfied. It is worth noticing, however, 
that when a small recoil is introduced the violation of the soft pion theorem is quite large 
in the full theory (see fig. 8). 



7 Conclusions 

In this paper we explored an approach that contains at its core the leading order in the 
heavy quark expansion and the next-to-leading order in ChPT, to derive expressions for 
the heavy-to-light pseudoscalar semileptonic form factors. The expressions are worked out 
in both standard and quenched ChPT. In the latter case we observe a familiar feature of 
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Figure 8: The ratio Fo(v-p)/(Jb/ f^) or f v /(f/f), which satisfies the soft pion theorem at zero recoil 
(v-p — > 0). The figure shows the ratio in the dynamical theory at which the violation of the soft pion 
theorem grows fast with the recoil momentum. The illustration is provided for three momenta and for 
masses corresponding to r <E (0, 0.5), i.e., m w e (0, 0.5) GeV. Notice that in the quenched theory this ratio 
does not depend on (v-p) and is equal to 1. 



quenched calculations, namely, the appearance of divergent chiral logarithms (quenched 
chiral logs), which makes it harder to compare the results with the formulas obtained in 
standard ChPT for B — > it form factors: one should find a fiducial window in which the 
masses are not so light so that the quenched chiral logs do not dominate the QChPT 
expressions, yet small enough for the chiral expansion to be meaningful. Furthermore, for 
the numerical estimates a number of low energy constants must be fixed by using both the 
available experimental data and the results of quenched lattice simulations. Our numerical 
estimates in which we use the QChPT expressions are highly sensitive to the value of the 
parameter g' (the coupling of the doublet of heavy-light mesons to a light rf meson) . Until 
a lattice computation of that parameter is made, we are not able to make firm quantitative 
statements. Even information about the sign of g' would be welcome. It turns out that, 
for g' = —g (which, on the basis of the large N c expansion, is the limiting value), we get 
a less pessimistic scenario: one can even find combinations of (v-p) and the light meson 
mass nip such that the quenching errors on F + ^(q 2 ) vanish. We were also able to find 
the combinations of (v-p) and nip such that the quenching errors do not depend on the 
value of g'; at these points and for small recoils the quenching errors are between 15% and 
25%. We reiterate that the numerical estimates do depend on the specific choice of the 
low energy constants. 

As for the form factor F (q 2 ), the present approach suggests that the quenching errors 
are large regardless of the value of the coupling g', the quenched value being larger than the 
unquenched one. Only at the point corresponding to zero recoil are those errors reasonably 
small. Away from that point, they are large. 
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The same observations apply also when the final meson is a kaon. In that case, by 
using the QChPT expressions, we were able to verify that the form factors obtained for the 
kaon consisting of degenerate and nondegenerate quarks are effectively indistinguishable, 
provided the value of the coupling is g' = —g. For g' = +g these uncertainties also become 
large. 

Finally, the formulas presented in this work may be useful in assessing the systematic 
uncertainty due to chiral extrapolations. We showed how to include the chiral logs to 
extrapolate from the region in which the pion is heavier than 400 MeV. As could have 
been anticipated, the estimated uncertainty of the chiral extrapolation depends on the 
mass from which the chiral logs are included in the extrapolation. 

We also provided the quenched chiral log coefficient, which might be useful if lattice 
calculations are performed with very light mesons. 

We verified that the ratio FQ^ 7r /(f B /f 7T ) satisfies the soft pion theorem, i.e., it is equal 
to 1 at zero recoil, in both theories. In the quenched theory that value remains unchanged 
even when a small recoil momentum is introduced. In the dynamical theory, instead, this 
ratio is significantly larger than 1 away from zero recoil. 
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Appendix A: Chiral loop integrals 

In this appendix we list the dimensionally regularized integrals encountered in the course 
of the calculation. For more details, see [52] and references therein: 



where 



ifj, 



(27r) 4 " e q 2 - m 2 16tt 2 



him) , 



1 1 



(q 2 - m 2 )(q-v - A) 16tt 2 A 



I 2 {m, A) , 



(40) 



m 



him) = m 2 ln^— -J — m 2 A , 
him, A) = -2A 2 ln(^) -4A 2 f(^) +2A 2 (1 + A) , 



(41) 



where A = 2/e — 7 + ln(47r) + 1. The function Fix) was calculated in ref. [53], for both 
negative and positive values of the argument: 



x 



VT^ 2 

x 

Vx^l 

X 



7T 



— tan 1 



x 



X/ 



In 



(x + Vx 2 



N < 1, 

\x\ > 1 . 



In addition to the integrals (40), one also needs the following two: 



■ . f d ""1 



q r 

(27r) 4 " e iq 2 - m 2 )iq-v - A) = IQ^ 2 

■4- e „ q n g u l 



iq 2 — m 2 )iq-v — A) lQir 2 



with 



him, A) = (-m 2 + ^A 2 ) In + ^(A 2 - m 2 )F 



;T2 j 



^A 2 (l + A) + im 2 (2 + 3A) + ^m 2 -lA 2 



Mm, A) = (2m 2 - ^A 2 ) In - 1(4A 2 - m 2 )F (£) 



+ ^A 2 (l + A) - ^m 2 (l + 3A) - ^m 2 + lA 2 



(42) 



[him, A) + him)] , 
A[him,A)r]^ + him,A)v^v v ] , (43) 



(44a) 



(44b) 



The functions Ji(m, A), him, A) differ from the ones in ref. [22] by the last two terms in 
(44) which are of (9(m 2 , A 2 ). These additional (finite) terms originate from the fact that 
is the (4 — e)-dimensional metric tensor. 
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Appendix B: Explicit expressions for the one- loop chi- 
ral corrections 

As already mentioned in the body of the paper, the chiral loop corrections to the form 
factors can be written in the form 



where the sum goes over all the graphs depicted in fig. 2. In what follows we give the explicit 
expressions for both form factors graph by graph and in both chiral theories (quenched 
and standard). 

B.l Quenched theory 

In the calculation of one-loop contributions, we made several approximations in order to 
simplify the final expressions. We make use of the fact that v-p > A* for the B — > 
P transition and thus consistently neglect the mass differences between B, B*, B S) and 
B* mesons in the loops. This neglect induces a spurious singularity in the expression 
for the diagrams (7a) and (7b) in fig. 2, at v-p — > 0. To handle those singularities we 
follow the proposal by Falk and Grinstein [25] and resum the corresponding diagrams and 
then subtract the term that would renormalize the B* meson mass. We recall that a (b) 
superscripts distinguish the diagrams without (with) hairpin insertion. The formulas for 
the Bj — > Pij transition (with the valence quark content of mesons Bj ~ bqj and Py ~ QiQj) 
are expressed in terms of the pseudoscalar meson mass Mf = 4fi mi: 




(45) 




v-p 3 



1 2tt 



5f (13a) 



18(4tt/) 2 



T 




+ 



a M; 
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(46) 



where T = for % = j, and T = 1 otherwise. The functions Ii(m), Ji(m,A) are given in 
appendix A. As for the form factor f v {v-p), the nonzero chiral loop corrections are 



T 



l2(Mj,vp) - h{M h vp) + - [Ii(Mj) - hiM,)} 



1 



6(4tt/) 2 (M 2 -M 2 



(aoMf-m 2 ) [I 1 (M j )+2I 2 (M j ,vp)} 



- (a M? - ml) (h(M t ) + 2I 2 (M i ,vp)) 
(a + (a M? - ™o)^f) (h(M t ) + 2I 2 (M h vp)) 



5f (Ua) 

6f (u b) 



2v / 6(4tt/) 2 
1 



(a M 2 - m 2 ) 



18(4tt/) 2 

+ (a Mf - ml) 



hiMj) + dh(M j ) 



Mj - M 2 



h{Mi) + 8h{M i 



M 2 - 



dM} 



a ih(M t ) + h(M 3 )] 



(47) 



In the wave function renormalization factors Z BP) we separate the one-loop chiral contri- 
bution SZ l g°p from the pieces coming from the counterterms SZp^p, i.e., 



Zb,p = 1 + 8Z B ,p = 1 + SZ l p°p + SZpfp 



(48) 



More specifically, 



r /7I00P _ ^ 

5Zb * ~Wff 



(2g*a M* - Ggg'M 2 - g 2 ml) In 0.) , 



+ a g 2 M 2 - mlg 2 



2g 2 M 2 a + 6gg'M 2 + g 2 m 2 ^j A (49a) 



5Z l ° op 



In (M 2 /M 2 ) 
9(4vr/) 2 l M 2 -M 2 



[2a M*M? - ml{M) + M?)] 



+ 2m2-a (Mi 2 + M|)} , 
while the counterterms contribute as 



5Z 



c.t. 



kim~ , 



5Z P ^-8L 5 ^-(m i + m j ) 



(49b) 



(50) 
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Thus the wave function renormalization factor for ir and K mesons reads 

P 

I f In (m s /m g ) 



1 - 8L 5 — 2m , 



Z K = 1 



9(47r/) 2 \ (m 2 K - m 2 



+ 2a m 2 K - 2mo| - SL 5 ^-(m q + m s ) , 



(51a) 



(51b) 



where we ignore the isospin-breaking effects and set m u = rrid = m q . 

Finally, we also display the expression for the heavy-light meson decay constant: 



a 



1 + 



6(4tt/) s 



a - ifV6Vl(0)} h{Mi) + (a M 2 - m\ 



dh(Mi 



dMf 



(52) 



in agreement with refs. [15, 16]. 



B.2 Full (unquenched) theory 

In this subsection we present the expressions for the form factors in the full theory. The 
nonanalytic contributions to the form factors in this theory have already been calculated 
in ref. [25]. Our results also include the analytic terms. As in the quenched case, we work 
in the isospin limit m u = m d = m q and neglect the differences of heavy meson masses in 
the loops. We now list the results for 8fv# in Bj — > mediated by the (V — A) operator, 
where P^ stands for the light pseudoscalar meson with the valence quark content q^qj. 
The form factor f p (v-p) receives the following one-loop corrections: 



6f^ = 



3 9 2 J \- r (7o) 



T < A 1 2ff 3 

Jx{m P i,v-p) 7T m P' 

v-p 3 



{ ^jf{ E C b?p>p, [Mmp^v-p] - ±f ml] } , 



1 2tt 
p 



1 



(Anf) 2 
1 



(53) 



where the coefficients CbjP'p^ depend on the final and initial states. A detailed list of 
coefficients includes the following. 



o FOR THE B —> K TRANSITION, 

(j(7a) _ n n( 7a ) — o r<( 7a ) 



'BttK 

r {l2a) _ _l ^(120) _ _1 r (l2a) 
BttK ~ 4' BKK ~ 2' ° B V K 



o, 



^ BKK — z ' 



/~<y' a l — O / Q. 
^BqK — 



«nr(9a) _ n 
^BttK — u > ^ BKK 



y(9o) 



0, 



1. 

'4' 



c 



(13a) 
BttK 



r (9a) 
^BriK 



n >nr(13a) _ _-, /-/(13a) 
u > BKK ~ L i Br/K 
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o FOR THE B — >• 7T TRANSITION, 



c 



(7a) 3 

Bn-rr 2 ' 

^(120) _ _2 

^ Btt-k 3 



= c 



(7a) 
(12a) 



1. 



c 



(7a) 1. 



(12a) 
Brjn 
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^Btttt — 2' 
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C 



0; C 



o FOR THE B s —> K TRANSITION, 
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2' 
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✓-,(120) 
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BlTTT 
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6' B s tvK — U 
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Brjn 
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BKir 
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2 ' °B?77r 
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'2' °B s ))if 
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1 ^(130) 
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The non- vanishing one-loop chiral corrections to f v (v-p) form factor are 



5/ (14a) = 



(4tt/) J 
1 



p' 



where the process-dependent coefficients D B .pr P .. are as follows. 

o FOR THE B K TRANSITION, 



n(4a) _ q D [.^, _ 2 r )l 4a; _ , j^waj _ 1/4 n 
BtxK ~ U > BKK ~ Z > ^PnK — -^PttE" — L /% ^ 



,(4o) 



,(4a 



(14a) 
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BKK 



o FOR THE B ^ TT TRANSITION, 

<l = 2, Dfl = 1, <t = 0; = -5/12, 



■1/2, <S = 



-1/3, <? 



(54) 



o FOR THE B s — >• K TRANSITION, 



D { i:i K = 3/2 , D ™ = 1 , d 



(4a) 



D 



(14a) 



l(4a) 
B sV K 



1/2; <i = -1/4 , 



-1/12. 
1/12. 
-1/2 , 



-1/12. 



The wave function renormalization factors Z for B mesons in the full theory are 

zg 2 



Z B u>d 
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1 - 



(47T/) 2 

3g 2 



3 1 

-him^) + h(m K ) + -h(m v ) 



+ k x m q + k 2 (m u + m d + m s ) , 



(4tt/) 2 

while for light mesons we have 



2I 1 (m K ) + -h(m v ) 



+ fciffl s + k 2 (m u + m d + m s ) , 



J K 
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(4tt/) s 
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4/^0 



4/iQ 
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(55) 



3(4^)2 [ 2I i( m ^) + J i(^x)] - 8L 5-^2~ 2m - I6L4— (m u + m d + m s ) . (56) 
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As in the previous subsection we close the list of results by showing also the corresponding 
formulas for the heavy-light decay constants. We have 

f Bs = _ Unf)2 [ J i( m ^) + ^ J i( m ^)] + xim s + H 2 (m u + m d + m s ) + ^Z Bs ^ , 



V™^ 1 (W) 

/*..< = ^={l -(^[^H + \h{m K ) + l/.K)] 



+ Kim q + x 2 (m„ + m d + m s ) + ^5Z Bud } , 

(57) 



in agreement with the results of refs. [15, 16]. 
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